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We study theoretically electronic Mach-Zehnder interferometers built from integer quantum Hall edge states, 
showing that the results of recent experiments can be understood in terms of multiparticle interference effects. 
These experiments probe the visibility of Aharonov-Bohm (AB) oscillations in differential conductance as an 
interferometer is driven out of equilibrium by an applied bias, finding a lobe pattern in visibility as a function 
of voltage. We calculate the dependence on voltage of the visibility and the phase of AB oscillations at zero 
temperature, taking into account long range interactions between electrons in the same edge for interferometers 
operating at a filling fraction v = 1. We obtain an exact solution via bosonization for models in which electrons 
interact only when they are inside the interferometer. This solution is non-perturbative in the tunneling probabil- 
ities at quantum point contacts. The results match observations in considerable detail provided the transparency 
of the incoming contact is close to one-half: the variation in visibility with bias voltage consists of a series of 
lobes of decreasing amplitude, and the phase of the AB-fringes is practically constant inside the lobes but jumps 
by 7r at the minima of the visibility. We discuss in addition the consequences of approximations made in other 
recent treatments of this problem. We also formulate perturbation theory in the interaction strength and use this 
to study the importance of interactions that are not internal to the interferometer. 

PACS numbers: 71.10.Pm, 73.23.-b, 73.43.-f, 42.25.Hz 



I. INTRODUCTION 

Recent experiments 110 on electronic Mach-Zehnder inter- 
ferometers (MZIs) constructed from integer quantum Hall 
edge states have attracted a great deal of attention. In these 
experiments Aharonov-Bohm (AB) oscillations are observed 
in the differential conductance of the interferometer. The most 
striking results concern behaviour at finite bias voltage. The 
visibility of AB oscillations shows a series of lobes as a func- 
tion of voltage, while their phase is independent of bias, ex- 
cept near visibility minima where it changes sharply by it. Our 
concern in this paper is with the theoretical understanding of 
these experiments. 

The observations are interesting from several perspectives. 
First, as was quickly appreciated^ it is plausible that the ef- 
fects arise from electron-electron interactions, because be- 
haviour of this kind does not occur in a single particle model. 
In addition, more seems to be required than a simple treat- 
ment in which inelastic scattering leads only to decoherence, 
since approaches of that kind cannot produce multiple side 
lobes in visibility of AB oscillations with increasing bias. It is 
remarkable that electron interactions should have the distinc- 
tive signatures found in this system, since integer quantum 
Hall edge states are usually modelled in the low energy limit 
as a chiral Fermi gas of independent particles. 11 The experi- 
ments therefore appear to reflect interaction physics that is not 
captured by the standard, universal description, but is robust 
enough to appear in many devices of varying designs. A sec- 
ond reason for interest stems from current efforts 12-20 to study 
interferometry in fractional quantum Hall states as a probe of 
fractional or non-abelian quasiparticle statistics. Against that 
background it is clearly important to understand unexpected 
interaction effects in much simpler, integer quantum Hall sys- 
tems. A third reason for interest is that the phenomenon seems 
to be an example of coherent many -body physics in a quantum 



system far from equilibrium. It invites comparison with other 
non-equilibrium quantum problems, from the Kondo effect 21 
to cold atomic gases. 22 

The design of an experimental device working as an MZI is 
shown in Fig.Q] It uses the edge states of a two-dimensional 
electron gas that is in an integer quantum Hall plateau. (Most 
experiments have been done at filling factor v=2, but broadly 
similar results have also been reported at v=l). The edge 
states serve as electron waveguides and are coupled at quan- 
tum point contacts (QPCs), which act as beam splitters. Cur- 
rent between, for example, source S 1 and drain D2 is mea- 
sured as a function of the voltage difference applied between 
sources SI and S2. Interference fringes are observed as oscil- 
lations in the differential conductance, either when the mag- 
netic flux density is varied by a small amount, or when a side 
gate is used to change the interferometer area or arm lengths. 
The visibility and phase of these oscillations vary with volt- 
age in the fashion already summarised. A physical scale is set 




FIG. 1: Schematic view of the electronic Mach-Zehnder interfer- 
ometer, which consists of a Hall bar with an island in it. A two- 
dimensional electron gas in a quantum Hall plateau occupies the 
shaded region. One edge state propagates along the lower edge of 
the Hall bar, from source SI to drain Dl, and a second edge state 
propagates around the island, from source S2 to drain D2. Tunnel- 
ing between these two edge states takes place at two quantum point 
contacts, with amplitudes t a and t^, at the points indicated by dashed 
lines. 
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by the bias voltage at which the first minimum in visibility oc- 
curs: the measured value, about in the first experiment, 2 
corresponds roughly to the chemical potential increase re- 
quired to add one electron to an edge state with length equal 
to that of the interferometer arms: for arm length d and edge 
velocity vf, this chemical potential increase is Tivf /d. 

Theoretical studies of coherence in electronic MZIs started 
before these experiments. Early work treated dephasing aris- 
ing from a variety of possible sources: interactions within 
the interferometer; 23 a fluctuating classical fieldri 2 ^ voltage 
probes;— or coupling to an external quantum bath. 27 None of 
these approaches generates the subsequently observed lobe 
pattern in the dependence of visibility on bias voltage. A 
further calculation^ based on a microscopic treatment of 
the effects of long-range interactions and using bosonization 
combined with a perturbative treatment of tunneling at the 
QPCs, shows that non-monotonic variations in visibility can 
appear for weakly coupled edge states, but without captur- 
ing the features found experimentally. By contrast, studies of 
models with additional structure, involving either a counter- 
propagating edge mode 29 or the pair of edge modes that arise 
at filling factor show that resonances can appear in 

that setting, which lead to lobes in visibility similar to those 
observed. These results are encouraging, but (as we discuss 
in Section lVirn i the models involved seem to us insufficiently 
generic to account for all experiments. More recently, approx- 
imate treatments of the effects of interactions at v=l when 
edges are strongly coupled by QPCs have generated some of 
the behaviour found experimentally^!^ There is good rea- 
son to think that these calculations identify some of the rel- 
evant physics, but the approximations used are non-standard 
and their domain of validity is unclear. 

In this paper we set out a detailed treatment of interaction 
effects in MZIs at filling factor v = 1. The approach is micro- 
scopic in the sense that it is based on the standard Hamiltonian 
for quantum Hall edge states, 11 and does not involve external 
noise. Our main results come from the exact solution of mod- 
els which have one simplifying feature: interactions that are 
restricted to the interior of the MZI; a short account of this part 
of our work has been presented previously. 34 We also present 
work in three further directions. One of these is an elemen- 
tary solution of the two-particle problem, which is a simplifi- 
cation of ideas from Ref. |3J. We believe that this calculation 
provides a useful illustration of the essential physics behind 
the phenomena we are concerned with, which is multiparti- 
cle interference. A second direction is a careful analysis of 
the approximations involved in Ref. 32. The third direction is 
the formulation of perturbation theory in interaction strength, 
which allows us to assess the importance of interactions that 
extend beyond the interior of the MZI. 

The organisation of rest of the paper is as follows. The 
two-particle problem is addressed in Section lU and the gen- 
eral microscopic description of the MZI is set out in Section 
Hm In Section|IV]we show how models with interactions only 
in the interior of the interferometer can be solved exactly. In 
Secton[V]we use this approach to study interferometers with 
various interaction potentials. We present a extended discus- 
sion of Ref. 32 in Section[Vl] and develop perturbation theory 



in interaction strength in Section [VTIl We summarise our con- 
clusions in Section [Villi Some technical details of the calcu- 
lations are given in appendices. 

II. TWO PARTICLE PROBLEM 

In this Section we set out a pedagogical treatment of the 
two-particle problem that illustrates how electron interactions 
affect the visibility of AB oscillations in an MZI. We con- 
sider an interferometer having both arms of the same length 
d and a propagation velocity vp for electrons. Denoting their 
separation by s, the flight time during which both are inside 
the interferometer isr= (d — s)/vp ■ We take the two elec- 
trons to interact with a potential energy U when both are in- 
side the MZI on the same edge, but not to interact otherwise. 
For simplicity, we consider first the case in which the mag- 
nitudes of the transmission amplitudes at the two QPCs are 
t a = h = l/v2, giving results for the general case later 

We solve the scattering problem for an initial state in which 
both particles are positioned on the upper channel before the 
first contact. We evaluate the probability for one or both par- 
ticles to exit the interferometer in the lower channel by sum- 
ming all quantum mechanical amplitudes that connect the ini- 
tial state to a given final state. We regard the expectation value 
for the total charge transferred from the upper channel to the 
lower channel as the analogue for the two-particle problem of 
the current in the many-body, steady state case. 

To establish some notation, consider in the first instance 
single particle scattering, initially for one QPC and then for 
an MZI. Amplitudes for the four scattering processes at one 
QPC are shown in Fig. [2] and the two possible paths through 
an MZI between an initial state in the upper channel and a 
final state in the lower channel are shown in Fig. [3] The am- 
plitudes associated with these two paths are given by products 
of the amplitudes arising at each QPC. Taking the total cur- 
rent to be proportional to the transition probability between 
the upper and lower channels, one obtains the standard result 
given in the caption to Fig. [3] with oscillations in the current 
as a function of the AB phase $. 

Now consider the two-particle problem with an initial state 
as described, in which both particles are on the upper channel. 
Paths to a final state with both particles in the lower chan- 
nel are shown in Fig. |4] and those to a final state with one 
particle in each channel in Fig. [5] (Because there is no dis- 
persion, particles cannot exchange positions in the scattering 
process. They may therefore be treated as if they were dis- 
tinguishable, and this is reflected in the figures by the use of 
different colours for the paths of each particle.) Without inter- 
actions the amplitude for a given pair of paths would simply 
be a product of contributions for each particle. Interactions 
contribute additional phase factors e ~ lUT / n when both parti- 
cles propagate on the same channel inside the interferometer. 
The average charge transferred between the upper and lower 
channels in the scattering process is 

I = 2\A 2 \ 2 +2|A!| 2 = 1 + cos(f7r/?i)cos$, (1) 

where both terms are multiplied by factors of two, since 
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FIG. 2: Possible paths and associated scattering amplitudes at a sin- 
gle QPC with transmission probability 1 /2. The phase a of the trans- 
mission amplitude is unimportant for a single QPC but contributes to 
the AB phase in an MZI. 
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FIG. 3: Possible paths and associated amplitudes for a single parti- 
cle to propagate through an MZI from an initial state on the upper 
channel to a final state on the lower channel. The phase difference 
$ = /3 — a is the AB-phase arising from enclosed flux. The com- 
bined amplitude A for transitions between these states is the sum 
of contributions from the two paths: A = — i e _l(Q+/3) < /2 cos($/2). 
The total current is proportional to \A\ 2 — i [1 + cos $]. 



A?, describes two-particle transmission, while single-particle 
transmission with amplitude A\ can occur for either parti- 
cle. AB oscillations are represented in this expression by 
the term in cos <&. Their strength is modulated by the factor 
cos(Ut/K). We can take the interaction strength U to play 
the same role in the two-particle problem as bias voltage in 
the many-body system, since increasing bias leads to reduced 
spatial separation between electrons entering the MZI above a 
filled Fermi sea, which in turn increases the interaction energy 
between these electrons. 

The phenomenon can be summarised by defining the visi- 
bility of AB oscillations. Let I max and I mm be the maximum 
and minimum values of / as $ varies. The visibility is 



From Eq. (Q} we have V = | cos(Ut/Ti)\, and hence a lobe 
pattern in V as a function of U. It is also evident from Eq. (Q]) 
that the phase of AB oscillations changes abruptly by tt at 
zeros of V. 

We use the same approach to calculate the current for the 
case of arbitrary tunneling amplitudes t ai t b at the QPCs. The 
result for the average charge transferred is 

I = 2[T a R b + R a T b + 2(T a T b R a R b ) 1 / 2 x 

x [l-4(T a i? a )sin 2 (f7r/^)] 1/2 cos$], (3) 

where T a<b = 1 — R a , b = t 2 b is the tunneling probability. 
The phase of AB oscillations is shifted by interactions, being 

$ = $+arccos{[l-4(T a i? a )sin 2 (C/r/ri)]^ 1/2 cos(Ur/7i)}. 

(4) 

The visibility is 

V = V x [1 - i(T a R a ) ^{Ut/Ti)] 1 ' 2 (5) 
where Vo is the single-particle value 

2[T a T b R a R b ] 1 ' 2 



Vn = 



T a R b + R a T b 



(6) 



An important difference between these results and the ones 
for the many body problem that we present in Section[V]is that 
here V does not decay at large U. In other respects, however, 
the two-body problem is illuminating. In particular, while the 
behaviour of V is not affected by the value of T b except for 
a multiplicative factor, if T a ^ 1/2 the zeros of the visibility 
turn into finite minima and the jumps in the phase of AB os- 
cillations become smooth rises. In the limit that transmission 
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FIG. 4: (Color online) The four possible paths and associated ampli- 
tudes for two particles to propagate through an MZI from an initial 
state with both particles on the upper channel to a final state with 
both particles on the lower channel. The combined amplitude is 

A 2 = -ie" i(a+/3) [l + e- iC/T/R cos$]. 
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at the first QPQ approaches T a = 1 or T a = 0, modulations of 
V with U disappear altogether. These features are also present 
in the many body problem. One consequence is that the lobe 
pattern cannot be obtained at leading order from a calculation 
that is perturbative in tunneling. 



III. MICROSCOPIC MODEL OF THE MZI 

Our model for the MZI is sketched in Fig. [6] The Hamilto- 
nian is 



n = n k 



Hi 



(7) 



It has three contributions: Hkin is the single particle term for 
an isolated edge; Hint represents electron-electron interac- 
tions; and ti t un — tif un + ti\ un describes tunneling at the 
QPCs labelled a and b. 

We consider initially edge channels of length L with pe- 
riodic boundary conditions, then take L — > oo. Allowed 
wavevectors are k = 2-KUk/L, with n% integer. Fermionic 
operators ct and Ck v , which create and annihilate an electron 
with momentum k on the edge rj, obey standard anticommuta- 
tion relations {c^c^ ' ,} = Sk p S„„r. In coordinate represen- 
tation the field operator ip„(x), which annihilates an electron 
at position x on the edge rj, is 



1 

7r 



E 



ik.i 



(8) 



With this notation 

Hkin = 



f L / 2 „ 

fiv F 2_. / i>+{x)d x 'ip v (x)dx. (9) 

77=1,2 J — 
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FIG. 5: (Color online) As in Fig. |4] but with one particle on the 
lower channel in the final state. The combined amplitude is A\ — 

-ie- ia e- iC/r/R sin$. 




FIG. 6: Schematic view of the MZI. Horizontal lines represent edge 
states with propagation direction indicated by arrows. These edge 
states are connected by two QPCs, shown as vertical dashed lines, 
with tunneling amplitudes t a and t^. The arm lengths between con- 
tacts are d\ and di, and the chemical potentials in the incident chan- 
nels are and fi2- The three different regions of the interferometer 
discussed in the text are labeled using Roman numerals. 



We represent interactions within each edge using the same 
symmetric potential U(x,x') and neglect interactions be- 
tween electrons in different edges. Introducing the density 
operator p„(x) = (x)^(x), we then have 



Hi 



L/2 



U(x,x')p v (x)pn(x')dxdx' . (10) 



?)=12 J-L/2 

Finally, taking the QPCs to be point-like, we write 



ti a tun = « B e fa #-(0)&(0) +h.c, 



(11) 
(12) 



Here v a and are tunneling strengths, from which the quan- 
tum amplitudes t a and tt, can be calculated. As in the previous 
section, the AB phase due to enclosed flux is <i> = j3 — a. 

This model can be solved exactly when interactions occur 
only between pairs of electrons that are both inside the inter- 
ferometer, in the region denoted II in Fig. [6] We present results 
in Section [V] calculated using three different choices for such 
internal interactions. The first of these is simply a charging 
energy 



U(x,x') = 



< x, x' < d 
otherwise . 



(13) 



An interaction of this kind is standard in the theory of quan- 
tum dots and for an MZI was treated approximately in Ref.l32l 
To test the robustness of behaviour to changes in the form of 
interaction, we also obtain results for two types of interaction 
potential that vary with electron separation inside the interfer- 
ometer, taking 



U{x,x') 



U(x-x') 




< x, x 1 < d 
otherwise 



and either an exponential dependence 

U(x - x 1 ) = ge- a \ x - x '\ 
or a Coulomb form 

U{x - x') 



sj{x- x 1 ) 2 + , 



(14) 



(15) 



(16) 
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IV. EXACT SOLUTION 

In this Section we give a full account of the solution out- 
lined previously in Ref . [34J. We study the interferometer at fi- 
nite bias voltage by computing the quantum-mechanical time 
evolution of an initial state in which the single-particle levels 
of fikin for the two edges are occupied up to different chemi- 
cal potentials, /ii and ^2. The interferometer reaches a steady 
state at long times, and we evaluate observables in this state. 
The challenge of the calculation comes from the difficulty 
of treating both tunneling and interactions non-perturbatively. 
Each part of the problem can be described by a quadratic 
Hamiltonian, but the appropriate variables are different in the 
two cases: fermionic for tunneling, and bosonic for interac- 
tions. A model which has only internal interactions can be 
solved exactly because in these circumstances the effects of 
tunneling at QPCs and of interactions can be handled sepa- 
rately. In the following we consider only edge states at filling 
factor v—\ and an initial state at zero temperature, but both 
these restrictions could be lifted within the approach. 

To treat time evolution we use the interaction representa- 
tion, taking as the free part of the Hamiltonian 

Hq = Hkin + Hint , (17) 

and as the 'interaction' part Htun- In this representation the 
time evolution of the fermion operators is given by 

j> v {x, t) = e lkot/n ^{x)e- lflot/h . (18) 

The wavefunction of the system, which we denote at t = by 
|Fs), evolves with the S -matrix 

£(i)=Texpj~jT n t un{r)dr\ , (19) 

where T denotes time-ordering. We distinguish operators in 
the Schrodinger and interaction representations by the ab- 
sence or presence of a time argument. In Section IIV Al we 
also use operators in the Heisenberg representation, and we 
indicate these with a subscript H. 

The presentation of the remainder of the calculation is or- 
ganised as follows. The observable we are concerned with is 
the current through the MZI, and we derive a convenient form 
for the corresponding operator in Section llV Al The simplifi- 
cations arising in a model with only internal interactions affect 
the calculation of the ^-matrix of Eq. (fT9l . which we describe 
in Section IIV Bl To find the time evolution of fermion oper- 
ators we use bosonization, as set out in Section IIV CI After 
bosonization Ho is quadratic and may either be treated using 
scattering theory (Section HVDb or diagonalised using a Bo- 
goliubov transformation (Section [IV El l. Inverting our trans- 
formations, we arrive in Section lTV Fl at an expression for cur- 
rent at long times, written in terms of fermion operators in the 
Schrodinger picture, and show that this expression is suitable 
for numerical evaluation. We give results for different choices 
of interaction potential in Section IV1 



A. Derivation of the current operator 

Although it is usual to write Htun as in Eqns. (fTTT i and 
(fT2l i. this is a shorthand since at finite tunneling strength the 
fermion field is discontinuous at QPCs. For that reason we 
regularise Htun by considering QPCs of finite extent w, tak- 
ing us — > at the end of calculations. Then, for example, 

Ht un nas tne f° rm 
1 f w 

^■tun = - / dx[v b e i ^+(d 1 +x)Md2+x)+h.c.} . (20) 
w Jo 

The current operator can be found in the standard way from 
the number operator N\ = J p\(x)dx for electrons on the 
upper edge, by evaluating its commutator with H. This gives 
for the current at QPC b 

h = ^ / v b e tf3 ipf(di + x)Md2 + x)dx (21) 

and a similar expression for I a , the current at QPC a. To cal- 
culate the position dependence of fermion fields within the 
QPCs we introduce operators in the Heisenberg representa- 
tion, with the time dependence 

A H (t) = jnt/K^-iiit/h (22) 

The equations of motion for the fermion operators at the con- 
tact b have in the interval < x < w the form 

(d t + v F d x )ip^ H (d 1 +x,t) = i^e~ lP ^ H (d 2 +x,t), 

(d t +v F d x )ip+ H (d 2 + x,t) = z^e l/3 ^+ ff (di +x,t). 

This system of equations has in the same interval the solution 

%H( d v +x,t) = ^2 M W Tp+, H (d v ,,t-x/v F ), (23) 

Tj' = l,2 

where M = exp(ixE) and the matrix E is 

v _ v b ( e- lf} \ 
^ ~ wv F \ e* ) ■ 

An explicit expression for the matrix M is 

M _( cos #(2) ie~ t/3 sm9(x) \ 
_ ^ ie^sme{x) cos9{x) ) ' 

with 6(x) = xvb/wv F h. Taking the limit w — > and writing 
e = +0, we arrive at a relation between fermion operators at 
d n + e, just after the tunneling contact b, and those at d v i — e, 
just before the contact: 

4>M d n + £ ^= E S S*^'H( d V - M), (24) 
with the scattering matrix 
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In this way the reflection r 0) & = cos8 ab and transmission 
ta.b = sm9 a ,b amplitudes at the contact a, b are expressed in 
terms of the angles 6 a j, — v a ^/hvF- 

Substitution of Eq. d23l into Eq. (f2TT > and reversion to the 
Schrodinger picture gives for the current operator If, the result 

h = ev F {tl[p 1H (di - e) - p 2H {d 2 - e)] 

+ t b r b [ie lf) ^+ H (di - e)^ 2H (d 2 - e) + h.c.]}, (26) 

in which we show explicitly that all operators are evaluated at 
a point infinitesimally before the contact. In the following we 
will omit e. An expression for the current operator at QPC a is 
obtained from Eq. ( f26b by substituting for di 2 and replacing 
9 b with 8 a . The total current operator is then I = I a + I b . 

When considering expectation values, denoted by (. . .) or 
by omitting hats, it is useful to separate the contribution from 
QPC b into two terms, I b = I t 



(i) 



r(2) 



with 



7 b (1) = eu F tg<ftO*i )-&(*)) 
li 2) = ev F t b r b [te lfl (G 12 ) +h.c], 



(27) 



where G\ 2 — ipf (dijiji^di). The term R is sensitive to the 
coherence between edges while 7^ is insensitive. Since there 

(2) 

is no coherence between channels before contact a, J c , = 
and the contribution to the current from this contact is 



iP=ev F t 2 a {p 1 (0)-p 2 (Q)). 



(28) 



(2) 

The term responsible for AB-fringes in the current is D ' and 
our general task is to calculate (G12). 

In experiment the differential conductance Q = e dl/d/ii 
(with p 2 fixed) is measured at finite bias voltage V — (p\ — 
1^2) /e. Q oscillates with $, having maximum and minimum 
values (? max and Smin- The AB fringe visibility is defined as 

?max ^min (29) 



B. Evaluation of the S -matrix 

We require the action of the S'-matrix, Eq. ( fT9l , on the 
initial state \Fs). This state is represented by a product of 
fermion creation operators acting on the vacuum and we need 
to find how S(t) transforms the fermion operators. Evaluation 
of S(t) is based on our restriction of interactions to the inte- 
rior of the MZI. Specifically, separating Ht U n into parts %f Un 
and nl un due to each QPC, we find (see AppendixlAl that 



i'^iun( t l)i'^tun( t 2)} - 



(30) 



for ti > t 2 - This leads to a factorization of the S-matrix into 
the product S(t) = S b (t)S a {t), where S a {t) is the S'-matrix 
calculated using H^ un , and S b (t) using H b un . A second com- 
mutator (see again AppendixlAl) 

[G 12 (t 1 ),n b tun (t 2 )] = o, 



also valid for t\ > t 2 , ensures that 

[S b (t)]+G 12 (t)S b (t) = G 12 (t), 



(31) 



so an explicit form for S b (t) is not required in the calculation. 
Since QPC a acts before the interacting region, it is easy to 
evaluate S a (t) ( Appendix IB1: we have 

for any t\, t 2 > and so may omit time ordering. The action 
of S a (t) on fermionic operators is a rotation in the space of 
channels and can be written as 



i> v {x) = [S a (t)]+^(x)S a (t). 
For < x < vptwe find the transformation 

with the rotation matrix given by 

S a = 



-it a e 



-it a e la 



(32) 



(33) 



(34) 



C. Bosonization 

To compute the time evolution of operators in the inter- 
action representation under Ho we use bosonization. 35 This 
gives us an exact correspondence between fermion and boson 
operators via the bosonization identity 

$„(s) = (2iray 1/2 F v e t2 ?^ x e- l ^ x \ (35) 
where bosonic fields are defined as 

k (*) = - £ i^/qL) 1 ' 2 (e^ x b qv + h.c.)e-" a / 2 (36) 

and a is an infinitesimal regulator, which does not enter the 
final results. The plasmon creation and annihilation operators 
(which have q > 0) obey bosonic commutation relations 



They can be expressed in terms of fermions as 



b+ =i(2,/«i) V2 E 



C h+qri C kr] 



(37) 



(38) 



k— — c 



The commutation relations for the fields (j> v (x) (omitting 
terms proportional to 1/L: see discussion in Ref.l35l) read 

[(/> v (x),d v 4>, r (y)] = -2niS (x - y)S w . 

The Klein factors F v , which change fermion number by one, 
satisfy the commutation relations 



{F ri ,F+}=2S w , [N V ,F V/ 



-S iF 
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with the standard expression for the particle number operator and 



N„ 



E 



C kr, C kr, ~ Yl ( 
k— — oo 



1^4,10) (39) 



in the edge r\. Here the vacuum state |0) satisfies 



e0, fc < 



c k „\0) =0, k > 0. 

The commutators [F v , bkri 1 ] and [7V^, bk v '] are zero. The elec- 
tron density operator is given by 



p v (x) = - — d x 4> n (x) + N„/L. 

Z7T 



(40) 



Since Hq does not couple channels, in the following we re- 
strict our attention to a single channel and omit channel labels 
until we reach Section HV Fl The kinetic energy Hkin f° r a 
single edge in bosonized form is 



TlVl 



L ' 2 dx 



2 .l-L/2 2tt 



^ It \\i i 27r ^ 



iV(N + 1) 



(41) 



Similarly, "H^t is quadratic, and given by 

#int = J/ [ U(x,x')p{x)p(x')dxdx'. (42) 

Using this form of the Hamiltonian, our objective is to ex- 
press the time-dependent boson field <p{d, t) in the interaction 
representation, in terms the boson operators b q and b q in the 
Schrodinger representation. We set out two approaches to this 
calculation. One is based on the formalism of scattering the- 
ory. We use this to treat interactions for which we can obtain 
simple expressions for plasmon scattering phase shifts. The 
other is based on a Bogoliubov transformation. We use it to 
study Coulomb interactions. 



D. Scattering approach 

The theory of plasmon scattering in spatially inhomoge- 
neous systems of quantum Hall edge channels has been stud- 
ied quite extensively. An early treatment of a Hall bar is given 
in R ef.pq and a recent application to an MZI is described in 
Ref. [29[ For the model we are concerned with, consider the 
equation of motion 

iW^{x,t) = [^{x > t),H \. (43) 

We separate 4>(x,t) = 4>^ [x,t) + 4>^ (x,t) into a part 
cM - 1 (x, t), proportional to N, and another part cM 1 ) (x, t), in- 
dependent of N. They obey 

N f d 

ihd t $ {0) (x,t) = -ihv F d x (t> {0) (x,t)+i— / U(x,y)dy 

L Jo 



2tt 



U(x,y)d y ^ (y,t)dy (44) 



iWt&W (x,t) = -ihv F d x 4> w (x,t) 



2tt 



U(x,y)d y $W dy (45) 



with initial conditions <p(°) (x, 0) = and (f>^ (x, 0) = <j>(x). 
Our aim is to find the Green function for Eq. d45l l. 

The basis functions for a mode expansion of $W (x,t) 
obey the time-independent Schrodinger equation 



u p fp(x) = -iv F d x f p (x) 



2-kTi 



U(x,y)d y f p (y) (46) 



and satisfy the orthonormality relation 

,L/2 

/ f P {x)d x f*(x)dx = ~2m5 p 

J-L/2 

The Green function can therefore be written as 



(47) 



e (48) 



K(x,y;t) = —J2fp(x)dyf;(y) 
p 

and we have 



fL/2 

(/> (1) (x,t)= / K(x,y;t)^(y)dy. (49) 

J-L/2 

Interactions within the MZI generate a frequency- 
dependent phase shift 5 P for plasmons, and the form of f p (x) 
on either side of the interaction region is (neglecting a correc- 
tion to the normalisation that vanishes as d/L — > 0) 



f P 0*0 = - 



I J ( 



iqx x <q 
e l(qx-5 q ) x> d. 



With periodic boundary conditions at finite L, the allowed val- 
ues of q are fixed by the condition f p (—L/2) — f p (L/2), and 
these determine the frequencies ui p = v F q. At long times 
and for large L, the quantity we require, K(d,y;t), can be 
expressed solely in terms of these phase shifts as 



l f°° 
K(d,y;t) = — J dp, 



i(p[d-y-v F t}-S v ) 



(50) 



From this we obtain 



4>^\x,t) = Y J {z q {x,t)b q +}i. C .), (51) 

in which the coefficients at long times have the form 

Zq {d, t) = - (2n/qL) 1/2 e Md-v F t)-i6 q (52) 

The long time limit of cf>^(d,t), which we write as (j>o(d), 
can also be expressed in terms of the phase shifts, as 



l (d) = 2tt— lim ^ 
L g^o q 



(53) 
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E. Diagonalisation by Bogoliubov transformation 

An alternative approach is to diagonalise the Hamiltonian. 
Substituting Eq. (f36t into Hq we obtain 



H = hv F Y q&tbq + 1/2) 

Y \fqk[u- q ,- k b q bt 



2L 



k,q>0 



Uq,-kbqK ~ u -q,kb q b k ] 



H N ^ + H N (54) 
where the matrix elements of the interaction potential are 



u q ,k = ^J o J e-^+^U (x, y) dxdy . 

They obey the relations u* k = and u q ^ k = U-k,-q- 

The last two terms in Eq. d54b involve the number operator 
N. The first of them appears because of our choice of nonuni- 
form interactions, which leads to a coupling 

H N4> = ih(N/L) (2n/L) 1/2 Y Vk[u Q , k b k - U ,- k b+] 

k>0 

(55) 

between the plasmon and the number operators: by contrast, 
in a system with translationally invariant interactions there 
would be no such coupling. The other term, Hjy, has the form 

H N = ^^N{N + 1) + ^ [ U(x,x')dxdx'. (56) 



L 2 



Since interactions in our model are limited to the finite region 
of length d, the second term in Eq. 061 ) gives a correction to, 
for example, the equation of motion of the Klein factor that is 
small in dj L and so vanishes in the thermodynamic limit. We 
therefore omit it in the following. 

Contributions to the Hamiltonian linear in TV are removed 
by making the shifts 



b q = b q + a q N , 
with coefficients a q given by 

1 



VFa q + — Y Vk]q[ u qkCt-k - u q _ k a* k ] 



(57) 



fc>0 



2tt 
L \~q~L 



1/2 



l q0 ■ 



(58) 

The Hamiltonian, Eq. d54l >. written in terms of these shifted 
operators, is diagonalised using a Bogoliubov transformation 
of the form 

p+ = J2( A Pi~K + B Pi~ b i)- < 59 > 

To preserve the commutation relations we require 

- B pk B+ q ) = S pq (60) 

fc>0 

J2(B pk Al q ~A pk Bl)=0, (61) 



which can be written in the matrix form 

A B \ ( A+ -B T \ i I 



B* A* 



-B+ A 1 



o i r 



(62) 



where / is the identity matrix. 
The result is 

H b = Y huj p(PpPp + V 2 ) + H N + const. (63) 

The time dependence of the transformed boson operators is 
given in the usual way in terms of their frequencies oj p as 
$ p (t) — e~ lUpt f3 p . Expressions for the coefficients A pq , B pq 
can be found from the commutator 

dt$ q = -i0 q , H b ] = -ifiwqhq (64) 
which leads to the linear system of Bogoliubov equations 

(ojp - v F q)A pq = \fqk[u qk A pk + u q - k B pk ] (65) 

k>0 

(lu p + v F q)B pq = - 2J \/qk[U- g>k A pk + U- q - k B pk ) . 

k>0 



From Eq. (1621 we can obtain the inverse of the Bogoliubov 
transformation, which we write in the interaction representa- 
tion as 

^(*) = ~ fc-^Bpq). (66) 

p>0 

Substituting d59l and (l57l i into equation ( l66*l l we obtain the 
time dependence of the bosonic fields in terms of the original 
operators b q written in the Schrodinger representation as 

4> (x, t) = $W (x, t) + 5^(* g (aj, t)b q + h.c), (67) 
where 

z q (x,t) = Y / (A; q f P (x)e-^ t 



S M /;(x)e^ 4 ) (68) 



p>0 



with 



f P ( X ) = ~Y ^h L ) (A pq e iqx - B pq e 

q>0 



%qx \ 



(69) 



k>0 



The term <f>(°> (x, t) in Eq. d67l >. arising from the operator 
shifts, is given by 

(o) {x, t) = 4> (x) - N Y^Zq. {x, t) a q + c.c.) (70) 

<;>0 

with 

Mx) = -NY(^/<lL) 1/2 {a q e^ -fee). 

q>0 

It is easy to check that at t = the field d67l i is equal to 
the Schrodinger operator (j>{x) of Eq. ([35]), since z q (x, 0) = 
- {2ir/qL) 1/2 e^and 0<W (x, 0) = 0. 

To make use of these results, the Bogoliubov coefficients 
A pq , B pq , a q and frequencies uj p are required. They can be 
found from Eqns. d65l ) and d58l , using a numerical treatment 
with a momentum cutoff. 
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F. Evaluation of the correlators 

In this section we explain how the treatment we have de- 
scribed of the bosonized Hamiltonian enables evaluation of 
the correlator G\2{t), which appears in Eq. ( f27l > and de- 
termines visibility of AB oscillations. When written using 
fermion fields in the interaction representation, it is 

G 12 {t) = (Fs\[S a (t)} + ^+(d 1 ,t)Md2,t)S a (t)\F S ) . (71) 

Our objective is to express ifj^(d\,t) and ^2(^2, t) in terms 
fermion fields ^ v {x) and if)i{x) in the Schrodinger repre- 
sentation, so that the expectation value in the state \Fs) can 
be computed. To this end we use the bosonization identity, 
Eq. (135) . and simplify notation by defining 



Then 



(72) 



(73) 



The bosonic field cj) v (d v ,t) which appears here in the inter- 
action representation is related at long times to one in the 
Schrodinger representation by Eqns. (l52l and (1531 . In addi- 
tion, the time evolution of the operator J- n (t) can be found in 
the usual way, via its commutator with Hq (the contribution 
from Hn4> is small in d/L and can be omitted in the thermo- 
dynamic limit), giving 



F v (t) = {2-Ka)- 1/2 F v e l ^ N ^ d ^ VFt) . 



(74) 



We wish to substitute for the Klein factor F v in this ex- 
pression. Consider a fermionic operator ifj v (z v ) with z n = 
d v — vpt. In the bosonized form, from Eq. ( 1351 ) it is 

= (2TTa)- l/2 F 11 e i ^ fl ^- VFt) e- l ^ ) . (75) 

Multiplying Eq. d75l > by from the right we obtain 

Trj (t) in terms of the fermion operators and bosonic fields in 
their Schrodinger representation, as 



(76) 



where the second equality holds due to the commutation rela- 
tions of Klein factors with bosonic operators. Substituting Eq. 
(|76*| | into Eq. (l73l we obtain 



ip v (d v ,t) = , 
which can be written as 



i> n {d ri ,t)=e- l Vr>e 



— p-i'Pvf>- i {<f>n( d v<' t )-'f>'i( z '))] 



4>n( z n) 



where the constant phase shift (p v is given by 

i f d<? -is 



(77) 



(78) 



(79) 



Finally, we substitute for b qv and 6+^ in 4>r)(x) in terms of 
fermion operators, using Eq. (1381 . with the result 



(80) 



Here the phase operator Q v is 

Qr,= / Q V (X - Z v )p rj (x)dx, 



(81) 



where kernel Q n (x) = J dqQ n (q)e lqx has Fourier transform 



Hi) 



--(e 

q 



1) 



(82) 



Eq. d80b is a key result which has a direct physical inter- 
pretation. An electron passing through the interferometer ac- 
cumulates a phase due to interactions with other electrons. 
This phase is a collective effect and it is represented at the 
point where the electron leaves the MZI by the operator Q v in 
Eq. (f80b . Contributions to the phase from the interactions with 
particles at a position x from the QPC b have a weight deter- 
mined by the kernel Q v (x). The form of the kernel is illus- 
trated in Fig.|7]for the case of a charging interaction, studied in 
Section IV Al The kernel has a maximum near x = 0, showing 
that interactions with nearby electrons are most important, but 
the phase is influenced by all the electrons which have passed 
the interferometer, although with contributions which decay 
with the distance x. The precise form of the kernel depends 
on the nature of the interaction potential and reflects the full 
many-body physics of the problem. A similar kernel appears 
in Eq. (1 1) of Ref. 32, but with a simpler form because of the 
approximations employed there. It is shown for comparison 
in Fig.|7j see discussion in Section lVTl 



o 




FIG. 7: The kernel Q(x) of Eq. d S 1 b for different values of the inter- 
action strength: 7 = 0.1 (short-dashed line), 7 = 0.2 (dot-dashed 
line), 7 = 0.5 (long-dashed line), 7 = 1.0 (dotted line) and 7 — 3.0 
(full line) compared to that of Eq. (1 1) of Ref. H3 at 7 = 1.0 (thin 
full line) 

Substituting Eq. (TSQb into Eq. (TTTT > we arrive at 

(G12M) = e^(Fs\[S a (t)]+^+(zi)e lA Mz2)S a (t) \Fs) . 
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Here R = Qi - Q 2 and $ = <p x ~ P2- The action of S a (t) + 
and S a (t) on the operators they enclose is given by Eq. d33l) . 
After this transformation the correlator reads 



1 

(G X2 {t))=j*- 5 i»e" 

k.q— — oo 



-ikz\+iqz2 . 



SM fl I dxQ 2 (x-z 2 )ip+(z 2 )M^)- (84) 



{F S \[c + ka e^c q p\Fs), (83) 

with summation over repeated indices a, (3. Here 1Z = 
[S a {t)] + RS a {t) is the rotated kernel 

il = S?*S?p I dxQ x (x - z 1 )^+(z 1 )^(z 1 ) 

Now evaluation of (Gi2(t)) reduces to the calculation of cor- 
relators of the form 

C m = (Fs\c+ expiiJ^MapctcMFs), (85) 

a/3 

where the indices specify both channel and momentum, and 
the matrix M is obtained from 1Z. One can show (see Ap- 
pendix that C m = D~ ^ det D with D constructed from 
the matrix elements of exp(iM) between the single-particle 
states that are occupied in the Slater determinant \Fs). We 
calculate C m numerically, achieving convergence of the re- 
sults when keeping up to 10 3 basis states and 400 particles in 
each channel: further details are given in Appendix iDl 

V. RESULTS FOR VARIOUS INTERACTION POTENTIALS 

In this Section we apply our theory to study interferome- 
ters with the three types of interaction potential introduced in 
Section |nil In the absence of interactions the interferometer 
at finite bias is specified by four dimensionless parameters: 
the tunneling probabilities t\ and t\ at the two QPCs; the di- 
mensionless bias eVy / did,2/2TrtiVF; and the ratio of the arm 
lengths d 2 jd\. The tunneling probability t\ at the second con- 
tact QPC affects only the overall scale for visibility of AB os- 
cillations, and we set it to t\ = 1/2. Interactions in general 
introduce another parameter, characterising their strength. Ex- 
ponential and Coulomb interactions also depend on a further 
parameter: the interaction range or the short-distance cutoff, 
respectively. 

A. Charging interaction 

Consider first the charging interaction, Eq. (I13t . It is 
characterized by the single dimensionless coupling constant 
7 = gd/2irhvp. Solving Eq. (l46l i we find 



/pW = " -7 




x < 
< x 
x > d 



with q = uj p /vf- Matching f p (x) at x = and x = 
d gives s p = (1 + tp)^ 1 and r p = t p s p , with t p — 
(g/2'KihiUp)(e iu >'' d/vF - 1). The phase shift 5 P of plasmons 
due to the interactions is 



^ — iS 



Similarly, we find 



» = (i + t;)/(i + t p ). 



4>o (x) = 2irjNx/L 



(86) 



(87) 



for < x < d, where 7 = 7(1 + 7) 1 . The contribution 
— ■£-d x (f>Q (x) — —jN/L to the density inside the interfer- 
ometer represents charge expulsion due to interactions: in the 
limit of strong interactions the average density inside the in- 
terferometer in the stationary regime is pinned at zero, inde- 
pendently of 7Y. 

The plasmon phase shift, as shown in Fig. [HJ varies lin- 
early with frequency at low frequency and falls to zero at 
high frequency. The maximum occurs at a frequency that in- 
creases with interaction strength and (for general interactions) 
depends on the shorter of two lengths: the interaction range 
and the arm length. For strong interactions, the phase shift at 
fixed frequency approaches the limiting value 5 p = Lu p d/vF- 
It then exactly cancels the kinetic phase qd. This remark- 
able cancellation together with the charge expulsion results 
in behaviour independent of arm length when interactions are 
strong. A similar cancellation was found for a different model 
in Ref.lM 




s p e iqx < x < d 



FIG. 8: Frequency dependence of the plasmon phase shift for a 
charging interaction with strength 7 = 0.1 (dot-dashed line), 7 = 
0.2 (long-dashed line), 7 = 0.5 (short-dashed line), 7 = 1.0 (dotted 
line), and 7 = 3.0 (full line). 

We now turn to behaviour of the interferometer. We con- 
sider first the case t\ = 1/2 and d x = d 2 - Results for the 
visibility of AB fringes as a function of bias voltage and inter- 
action strength are presented in Fig. [9] Without interactions, 
visibility is independent of bias, having a value fixed by the 
tunneling probabilities in the QPCs. Interactions generate a 
dependence of visibility on bias. Visibility at small bias is un- 
affected by interactions, because in this regime there is only 
one extra electron inside the interferometer at a time, but with 
increasing voltage the visibility follows a sequence of lobes 
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1.5 




1 2 3 4 5 6 

eVd/27r/lVF 



FIG. 9: Visibility as a function of bias voltage for an MZI with charg- 
ing interactions, and with di — d^ and 1? a = t% — 1/2, at interaction 
strengths: 7 = 0.1 (short-dashed line), 7 = 0.2 (dot-dashed line), 
7 = 0.5 (long-dashed line) and 7 = 1.0 (full line). The phase of the 
AB-fringes (not shown) jumps by n at zeros of the visibility. 



separated by zeros. The width in bias voltage of lobes is in- 
versely proportional to 7 for small 7 and saturates at a value 
close to irhvF/ed for 7> 1. The phase of AB-fringes is in- 
dependent of bias inside the lobes and jumps by ir at the zeros 
of visibility. Both features, the lobes in the visibility and the 
phase slips, match those of experiment by Neder et. al. (see 
Figs. 2 and 3 ofRef.0). 

It is interesting to note that the visibility can become larger 
than one. This signals a negative value of the differential con- 
ductance. Similar behavior has been reported previously. 28 

Results for an interferometer with unequal arm lengths, 
di/di — 1-2, are presented in Fig. [TOj The gross behav- 
ior at intermediate and large interaction strengths is similar 
to that for an interferometer with equal length arms. The 
visibility minima for d\ ^ d%, however, are not exact ze- 
ros: they approach zero at large interaction strength, but dis- 
appear altogether in the opposite limit of weak interactions. 
There is a corresponding evolution with interaction strength 
in the dependence of the phase of AB fringes on bias voltage. 
In the absence of interactions this phase varies linearly with 
bias for an MZI with different length arms because the Fermi 
wavevector kjr is linear in bias and the phase difference be- 
tween particles traversing the two arms is fci?(c?2 — di). With 
increasing interaction strength the phase dependence on bias 
develops into a series of smooth steps, each of height tt. The 
risers of these steps coincide with minima of the visibility. 
Strikingly, for the interacting system phase steps at minima 
of the visibility persist for d\ — g?2, even though in this case 
phase would be independent of bias without interactions. Be- 
tween risers the AB phase is almost independent of the ratio 
^2 / di when interactions are strong, because of the cancella- 
tion between kinetic and interaction contributions to plasmon 
phase, as discussed above. The stepwise phase variation we 
find at large interaction strength also matches experimental 
observations (see Fig. 2 of Ref.|j). 

Behavior of the visibility is insensitive to the transmission 
probability t\ at QPC b, apart from the overall scale. De- 
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FIG. 10: Visibility (upper panel) and AB phase (lower panel) as a 
function of bias voltage for an MZI with dz/di = 1.2 and t„ = tjj = 
1/2 at interaction strengths: 7 = 0.1 (short-dashed line), 7 = 0.2 
(dot-dashed line), 7 = 0.5 (long-dashed line) and 7 = 1.0 (full 
line). The linear dependence of the AB -phase in the noninteracting 
case 7 = is shown on the lower panel (thin full line). 



partures from t 2 a = 1/2, however, like unequal arm lengths, 
eliminate the exact zeros in visibility, leaving only sharp min- 
ima. The dependence of visibility on voltage for t\ — 0.75 is 
shown in Fig. [TT] 



B. Exponential interaction 

We next consider the interaction potential of Eq. ( I"f"5l >. 
which decays exponentially with separation. It is charac- 
terised by g, the interaction strength, and 1/a, the range. In 
this case the Schrodinger equation, Eq. d46l i. for < x < d 
has the form 

- iu p f p (x) + v F d x f p (x) = / e ~ aVx ^ d yfp(y) d V ■ 

(88) 

Differentiating Eq. d88l twice with respect to x we obtain 

d X xxfp(x) - i(uJp/v F )(d xx f p {x) - a 2 f p (x)) 

- (a/d)(27 + ad)d x f p {x) = (89) 
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FIG. 11: Visibility (upper panel) and the AB-phase (lower panel) as 
a function of bias voltage for an MZI with c?2 = di, t% = 0.75 and 
tl = 1/2, and charging interactions of strength: 7 = 0.1 (short- 
dashed line), 7 = 0.5 (dot-dashed line), 7 = 1.0 (long-dashed line), 
7 = 2.0 (dotted line) and 7 = 3.0 (full line). 



with 7 = gdj (2ttTivf) as before. This equation has a solution 
of the form 



Solving Eqns. ( l90t and d92l numerically we obtain the phase 
shifts from 



S p = -Axg[/ P (d> 



^—iud/vp 1 



(93) 
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FIG. 12: Frequency dependence of the plasmon phase shift for expo- 
nential interactions of strength 7 = 0.5 and range: ad = 2 (dotted 
line), ad = 1 (dot-dashed line), ad — 0.5 (dashed line), and for a 
charging interaction, ad — (solid line). 



The frequency dependence of this phase shift is shown in 
Fig. Q2] for different values of the interaction range 1/a and 
the interaction strength 7=1. The main features, of a linear 
variation at low frequency and a phase shift approaching zero 
at high frequency, are independent of range. 

The resulting fringe visibility in an MZI with these inter- 
actions is illustrated in Fig. [T3] At fixed interaction strength 
the visibility has zeros at values of the bias voltage which are 
set by the energy scale TiVFOt for interactions with range much 
shorter than the arm length and by 2irhvF/d for interactions 
with range in the opposite limit. 



fp(x) — A r , 



(90) 



C. Coulomb interaction 



where the A n are in general complex coefficients and the 
wavevectors k n are obtained by solving the cubic equation 

u) p = kv F [l + U(k)/2iThv F ]. (91) 

Here U(k) = 2ag/(k 2 + a 2 ) is the Fourier transform of the 
potential in Eq. lfI5V Substitution of Eq. (J90J into Eq. (l88l l 
gives two linear equations on A n and the boundary condition 
/p(0) = 1 yields a third equation. These determine the coef- 
ficients A n and read 

n— 1 n—1 

3 , 

and y^A n - — — =0. (92) 



Finally we consider the unscreened Coulomb interaction, 
Eq. ( [ToT l. which is characterized by its strength j c = 
g c /2irTiVF and the short-distance cutoff a c . We treat the 
regime a c <C d; in the opposite limit it is similar to the charg- 
ing interaction discussed in Section IV Al (Coulomb interac- 
tions have been studied previously, 28 without our restriction 
that they act just within the MZI, but only perturbatively in 
tunneling at QPCs.) To calculate the plasmon phase shifts in 
this case we solve the Bogoliubov equations [Eq. d65l ll nu- 
merically. Results for the bias dependence of visibility are 
presented in Fig. [14] The visibility again shows lobes. Their 
width in bias voltage is inversely proportional to the interac- 
tion strength j c for weak interactions and is set by the inter- 
ferometer energy scale 2Trhv F / d for strong interactions. 

In summary, while the detailed shape of oscillations in visi- 
bility with bias voltage depend on the model used for interac- 
tions, the main features are independent of this choice. 
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FIG. 13: Visibility as a function of bias voltage for an MZI with 
exponential interactions [Eq. 115H . with d,2 = di, — t'l = 1/2 
and 7 = 0.5, for interaction ranges: ad — 0.5 (dashed line), ad = 1 
(dotted line), ad = 2 (solid line). 
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FIG. 14: Visibility as a function of bias voltage for an MZI with un- 
screened Coulomb interactions [Eq. J16H . taking a short-distance 
cutoff a c — O.ldi, c?2 = di, t„ = = 1/2, and interaction 
strengths: 7 C = 0.05 (dot-dashed line), j a = 0.1 (dashed line) and 
7 C = 0.15 (solid line). 



VI. DISCUSSION OF THE APPROACH OF NEDER AND 
GINOSSAR 

In a recent paper 32 Neder and Ginossar have presented an 
approximate treatment of interaction effects in an MZI, using 
the charging interaction of Eq. [13] They arrive at a result sim- 
ilar in form to the exact one given above in Eqns. (f8C)b and 
dHTJ, but with a different kernel Q(x): see the comparison 
in Fig. [7] Making further simplifications in the course of a 
numerical evaluation, they obtain nodes in the dependence of 
visibility on bias voltage. We discuss their approach in this 
section (we note that various choices of kernel and of approx- 
imations schemes in this type of numerical evaluation have 
also been discussed in Ref. 31). We show that their central ap- 
proximation is equivalent to neglect of the chiral anomaly in 
the Tomonaga-Luttinger model, and explain physically what 
this entails. We also examine how far the differences between 
results for visibility in Ref. |32] and those in the present pa- 



per are due to this central approximation, and how far they 
stem from simplifications of the numerical evaluation made in 
Ref. To do this we make a numerically exact calculation 
of the fringe visibility using the kernel of Ref. 32 and compare 
the outcome with our own results. 



A. Single edge with translationallv invariant interactions 

The approach under discussion starts from the Heisenberg 
equation of motion for an operator An(t). Consider first a 
single edge channel without QPCs. In this case the equation 
of motion is 

ihd t A H {t) = [A H (t),H ] (94) 

with initial condition Ah(0) = A, where A is the operator in 
the Schrodinger representation. In the absence of interactions 
Eq. (l94l for the field operator has the solution 

^+{d,t)=^+(e,t-d n /v F ). (95) 

Next, include translationally invariant interactions U (x, x') = 
U(x — x'), with U(0) = to avoid self-interactions. The 
equation of motion is 

ih(d t + vpd^^x.t) = 

/CO 
U(x- x')p H (x',t)dx'. (96) 
-OO 

It apparently has the solution 

V3+ (x, t) = e a *<?>H%{x ~ v F t, 0) (97) 
where the phase operator 9h(x, t) is given by 
t f°° 

9 H (x,i) = - U(x')p(x-v F t + x',0)dx'. (98) 

" J — OO 

This calculation is essentially a version in the Heisenberg 
picture of the solution of the Tomonaga-Luttinger model by 
Luttinger, 38 who used a canonical transformation to diago- 
nalise the Hamiltonian. It is exact provided the number of 
particles in the system is finite. The difficulty, of course, is 
that for finite particle number a model with linear dispersion 
has no ground state, and so one wants to introduce a filled 
Fermi sea. Unfortunately, as shown by Leib and Mattis^ the 
calculation is then no longer exact, because density operators 
which commute when particle number is finite no longer do 
so in the presence of the Fermi sea. 

To illustrate the difference in physical behaviour between a 
system with a finite number of particles and one with a Fermi 
sea, it is useful to consider the time evolution of an initial state 
in which single particle orbitals are all occupied for wavevec- 
tors k in the range < k < fcp, and others are all empty. 
We denote this state by \k F ) and compare its evolution with 
that of the state with a Fermi sea, in which all orbitals with 
k < k F are occupied. The latter is an exact eigenstate of H 
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for arbitrary choice of interaction, and so has trivial time evo- 
lution. By contrast, \kp) is not an eigenstate in the presence 
of interactions. To be specific, we calculate the time depen- 
dence of the particle number h(Q) = CqCq in an orbital with 
wavevector Q. Within a short time expansion we have 



(fi(Q,t)) = (n(Q))+it([H ,h Q }) - -{[H ,[H ,n Q ]])... 

(99) 

For the state |fcp) the average ([Hq, n(Q)}) = 0. The second 
order term is conveniently expressed in terms of 



where 0h(x, t) satisfies the equation 



U(q) = / dxU(x)e iqx 



(100) 



the Fourier transform of the interaction potential. Taking fc F 1 
to be much smaller than the interaction range, and writing 
Q = k F + P with P > 0, we find 

(n(Q,t)) = t 2 (H Q h(Q)H Q ) + ... 

t 2 r°° 

= — J p dqq\U(q)\ 2 + 0(t 3 ). 

Similarly, for — fcj? < P < we obtain 

t 2 f°° 

(n(Q,t)) = l-_ J ^dqq\U(q)\ 2 +0(t 3 ). (101) 

This shows that interactions scatter particles from the orbitals 
they occupy initially in \kp) to others of larger and smaller 
energy. The pair of particles involved in a typical scattering 
event has one initial wavevector just less than kp and the other 
just larger than zero. One particle is scattered to a state with 
wavevector k > kp and the other to a state with negative 
wavevector. The rate for this process remains finite even when 
kp is large, but in the presence of a Fermi sea these scattering 
processes are blocked by Pauli exclusion. 

B. Interactions confined to a finite region 

To demonstrate that no other approximations are involved 
in the derivation of the kernel of Ref. |32|, we next set out a 
calculation equivalent to the one leading to Eq. d97j >, but for 
interactions that operate only in the finite region < x < d 
and are translationally invariant within this region, so that 



Hint — 



d pd 



Jo 



U{x — x')p(x)p(x')dxdx' 



(102) 



Then the equation of motion for ip + (x,t) can be solved by 
integrating forwards in time separately in each of three regions 
and matching at boundaries. For x < we have 

-4> + {x,t) = -4> + {x- vpt,0). (103) 
Within the approximation under discussion, for < x < d 

TP+(x,t) = e i8 »^i>+(x-VFt,0), (104) 



(dt + v F d x )9 H (x,t) = - J U(x-x')p H (x',t)dx' (105) 
with the boundary condition 8(0, t) = 0. This has the solution 



H (x,t) ^^—i (x' + v F t)+ 

nv F U_ VFt 

■'■ / + / (d + x) } x 



x—vf^ J d—Vpt 



U(x)p H (x',0)dx', (106) 



in which x, = x — vpt — x'. 

In particlular, for x — d we have 



1 



H (d, *) = it- / x'U{d - x')p H (x', t)dy 



i r 2d 

+ - — / (2d- x')U(d- x')p H (x',t)dx' . (107) 

TiVF Jd 

In the region x > d we can use results from < x < d as a 
boundary condition to obtain 

i>+(x,t) =e i§H{d ' t) i)^(x-VFt,0). (108) 

For the charging interaction, Eq. ( fT3b . this gives an expression 
for the kernel Q(x) appearing in Eq. ( T8Tb which is equivalent 
toEq. (11) of Ref. IS 



C. Interferometer with internal interactions 

For an interferometer with only internal interactions, the 
foregoing discussion can be combined with the approach de- 
scribed in Sections TlV GUV Dl and llV El In particular, making 
use of Eq. ( I108l l in place of Eq. (1821) . one arrives at alternative 
results for visibility. We have evaluated these without further 
approximation for charging interactions. They are displayed 
in Fig. [15] In this approximation the visibility as a function of 
bias shows a lobe pattern similar to that given by our exact so- 
lution (Fig. [9]) but some significant differences are apparent. In 
particular, within the approximate treatment visibility at small 
bias is reduced from its value in the non-interacting system, 
and more so as interaction strength is increased. This loss of 
coherence mirrors the scattering described by Eq. ( 11011 ). It is 
not a feature of the exact treatment for a system with a Fermi 
sea: in that case at small bias extra electrons are dilute and 
pass independently through the interferometer. 



VII. PERTURBATION THEORY IN INTERACTIONS 

The approach we have used to obtain exact results can be 
applied only to models in which interactions act solely within 
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FIG. 15: Visibility as a function of bias voltage for an MZI with 
charging interactions, calculated using the kernel of Eq. (1 1) of Ref. 
l32l . with di — d,2 = 1 and t\ = t% — 1/2, for interaction strengths: 
7 = 0.1 (dashed line), 7 = 0.2 (dot-dashed line) and 7 = 0.5 (solid 
line). 



the MZI. In order to gain some understanding of the conse- 
quences of more general interactions, it is useful to formu- 
late perturbation theory in powers of interaction strength. We 
outline such a calculation in this section. Since details are 
quite messy (around forty separate terms appear at first or- 
der in interaction strength), we limit ourselves to sketching 
how a new physical effect - a second harmonic in the varia- 
tion of differential conductance with AB flux - appears when 
there are interactions between the interior and exterior of the 
MZI. This expansion in powers of interaction strength, which 
can be applied for arbitrary values of the tunneling amplitude 
at QPCs, is complementary to the perturbation theory in tun- 
neling strength, developed for arbitrary interaction strength in 
Ref.ii 

We separate the Hamiltonian, Eq. (|7), into the single- 
particle contribution Hi = Hkin + T^tun and the interaction 
term Hint, an d start from a single-particle basis of scattering 
states of Hi, These states are labelled by energy frvpq and by 
the channel r\ (denoted 1 and 2 in Fig. [6) from which particles 
are incident on the MZI. They have amplitude e lqx ip qv (x, n) 
at point x in channel k. Thus, for example, <p q i(x, k) takes 
the form given in Table 1 . 



K 


x < 


< x < d K 


d K < x 


1 


1 




r a r b -t a t b e i ^- d ^ +ts - a 1 


2 





~it a e- la 


-it a r b e~ la -irahe^^-^-V 



TABLE I: Amplitudes of the scattering state tp q \ (x, k). 

We consider an initial many -particle wavefunction |^o) m 
which scattering states are occupied up to Fermi wavevectors 
pi and P2 respectively in the two incident channels, represent- 
ing a bias voltage V = Tivp(pi — p<i)/e. Using the interac- 
tion representation A(t) = e lHlt ^ n Ae~ lfflt ^ n based on Hi 
(different, of course, from that employed in Section ITVT l. we 
evolve the initial state forward in time from t = —00 to t = 0. 
The current at zeroth and first order in ThLi n t is 



(/(t)) t=o = (*o|/(0)|fo} 

ft) 



Hint (t) , / (0) 



|*o) dr. (109) 



To evaluate the interaction term in this expression we take 
matrix elements in the basis of single-particle scattering states, 
denoting for brevity the pair of labels q a , n a by a. For an 
interaction U(x, x') that does not couple channels we write 



Uabcd 2 



dxdx' U(x, x') e *Ka«-9.)»+(ge-ft)»'] 

^l a K{ X )^*q b K{ X ')^on{x')ip qdK (x) (110) 



and define the antisymmetrised combination U a bbd — 
\{U a bbd — Uabdb)- We evaluate matrix elements I a t, of the 
current operator by considering in the first instance a tunnel- 
ing contact of finite width w, then taking the limit w —> 0, as 
in Section lTV Al For example, 

hi ,1 - evF (f a + r^^-*) _ t 2 a t 2 b e id ^ 



+t a r a tbn 



e i(a-/3+qd 1 -kd 2 ) _|_ e i(/3-a+<jd 2 -fedi) ~j 



With this notation, and denoting by n a the average occupa- 
tion of the state a, the first order term from Eq. (1109b has the 
form 



2n 3 h 



dr / dk a 



dk n 



x UabbJcaUb (n a - n c ) e »Hfc»-fcc)T , 

Evaluating the integrals on r, k a , kb and k c we obtain an ex- 
pression in which only integration on x and x' from Eq. dl 101 ) 
remains. Because the factors Lp qrj (x, k) take different forms in 
each of the regions I, II and III, defined in Fig.|6l these final in- 
tegrals naturally separate into distinct contributions according 
to the possible locations of each of two interacting particles. 

As an illustration of this general approach we consider the 
contribution to current arising at first order from the inter- 
action between a particle in region II and one in region III. 
This is made up of an exchange term Sip 3 and a Hartree term 
Sljf. To write expressions for these in a concise form we limit 
ourselves to the case di = d,2 = d and define the constant 
Tab = tahr a r b . Then 



5lf 



^- sin 2$ 



: / dx 
'0 



'dx>Uix,x>)^-^ X - ^ 



(x — x') 2 



and 



SI 23 



H 



leT 1 

— ^(pi-p 2 ) 2 sin2$ 
irn 



x / dx 
lo 



dx' U (x, x) . 
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As advertised, both terms involve the second harmonic sin 2$ 
of the AB phase They display the antisymmetry expected 
for the current under reversal of both V and $. A second har- 
monic is also produced by the interactions between the elec- 
trons in the regions I and III which we do not present here. 
By contrast, a similar perturbative treatment of interactions 
between electrons in region II generates only zeroth and first 
harmonics of $, as expected from the results presented in Sec- 
tion |IV] 



VIII. DISCUSSION 

In summary, we have argued that electron interactions are 
the origin of the experimentally observed dependence of vis- 
ibility of AB oscillations on bias voltage in electronic MZIs. 
This is illustrated most simply by calculations of two-particle 
interference effects, described in Section[n] and demonstrated 
in detail by exact results for the many body system, presented 
in SectionjV] Our calculations rely on a simplified form for in- 
teractions, but we believe our choice is quite reasonable. Our 
central approximation is to neglect interactions between an 
electron inside the MZI and one outside. In practice, such 
interactions will anyway be screened by the metal gates that 
define the QPCs. We also neglect interactions between a pair 
of electrons that are both outside the MZI. This is unimpor- 
tant: before electrons reach the MZI, such interactions do not 
cause scattering because of Fermi statistics, while after elec- 
trons pass through the MZI, these interactions cannot affect 
the current. We have given results for three forms of inter- 
action within the MZI, finding their main features to be in- 
dependent of details of the model. These features include a 
series of lobes of decreasing amplitude in AB fringe visibil- 
ity as a function of bias, with jumps by it in the phase of AB 
oscillations near minima in visibility. 

The width in bias voltage of the central visibility lobe de- 
fines an energy scale. For our model of charging interactions 
this scale is of order g at large 7. Taking vp = 2.5 x 10 4 ms _1 , 
d = 10/im and the permittivity e = 12.5 of GaAs, we estimate 
from the capacitance of an edge channel g ~ 10/ieV. This is 
similar to the experimentally observed value of about 14 /ieV, 
given in Ref . 

We close by comparing the results we have presented with 
those from other approaches. As a first step, we note that, 
while the existence of multiple side lobes in the visibility of 
AB oscillations as a function of bias voltage cannot be ac- 
counted for by a simple treatment of dephasing, a single side 
lobe can emerge from a simple, phenomenological treatment. 
To see this, consider the dependence of current I(V, on 
bias V and AB flux $. Assuming that there is a just one har- 
monic in AB oscillations, we have 

7(V,$) = I o (V)+h(V)cos($ + 0(V)). (Ill) 

With this notation, the visibility is 

[{d v h{v)f + {h{v)d v e{v)f]^ 

\dvh(v)\ ■ K ] 



It is reasonable to expect quite generally that h(V) increases 
with V for small V, has a single maximum, and decreases 
towards zero at large V, so that dy h (V) has a single zero for 
< V < 00. If in addition the phase 9(V) of AB oscillations 
is independent of V (and only in this case), the existence of a 
single side lobe in visibility follows. 

The theoretical difficulty, then, is to understand the obser- 
vation of multiple side lobes. Calculations to date that gener- 
ate such behaviour can be divided into three categories. One 
of these 29 involves a plasmon resonance between one arm of 
the MZI and a counter-propagating edge state at the bound- 
ary of the Hall bar. Such a coupling of the MZI to another 
edge state is not an essential part of the interferometer design, 
and in this sense the mechanism is not an intrinsic one. For 
that reason, it seems unlikely to provide the explanation for 
observations in many different samples of varying designs. 
The most important comparison is therefore between a sec- 
ond category of explanation, 30 which is based on coupling be- 
tween the two channels existing at each edge for filling factor 
v = 2, and the third category, which is formed by the calcu- 
lations that we have presented, together with earlier approxi- 
mate discussions 32 ' 33 of similar physics, and has been worked 
out for a system at v = 1 . 

According to the approach of Ref. 30, multiple side lobes 
should be found only at v = 2, which seems indeed to be 
the case experimentally. Some important discrepancies be- 
tween this theory and experiment remain, however. One is 
that, within the approach of Ref. [3(1 the shape of the enve- 
lope of the lobe pattern is controlled by the difference in in- 
terferometer arm lengths: in particular, for an interferometer 
with arms of equal length, visibility does not fall to zero at 
large bias. This is in conflict with observations. It requires 
one to assume 30 that two separate physical processes are in- 
volved: the process included in the theory, which leads to 
multiple zeros in visibility at v = 2, and another one, omit- 
ted from the theory (such as dispersion of the edge modes), 
which controls decay of the envelope. Moreover, if this the- 
ory is adopted at v = 2, the existence of a single side lobe 
at v = 1 must be attributed to a separate dephasing mech- 
anism, following Eq. ( II 121 i. In experiment, there appears to 
be a common voltage scale determining all aspects: the ze- 
ros in visibility at v = 2, the envelope of the lobe pattern at 
either filling factor, and the position of the visibility zero at 
v = 1. It would be a surprising coincidence if two separate 
mechanisms were both to involve the same scale. An expla- 
nation of multiple side lobes in visibility based exclusively 
on dispersionless slow and fast edge modes at v = 2 there- 
fore seems problematic. By contrast, the calculations we have 
presented generate zeros of visibility and an overall decaying 
envelope for visibility from a single mechanism. Depending 
on interaction strength we find either multiple side lobes or 
only a single prominent side lobe. The mechanism is essen- 
tially plasmon dispersion. It is likely that a full understanding 
of experiments will require a combination of both aspects - 
dispersion and the existence of two modes at v = 2. It is 
however, an important point of principle, demonstrated by the 
calculations we have described, that multiple side lobes are 
not an exclusive consequence of coupling between a slow and 
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a fast mode at v = 2. 

A definite qualitative feature of our results is that exact ze- 
ros of fringe visibility at certain values of bias voltage appear 
only for an interferometer with equal length arms, in which 
the transmission probability at the first QPC is tuned precisely 
to the value one half. Changes in transmission probability and 
(though to a much lesser extent, if interactions are strong) arm 
length from these values convert exact zeros to finite minima 
in visibility. Such sensitivity to transmission probability has 
not been reported experimentally. The reason it appears in our 
theory can be understood starting from Eq. (II 12b : exact zeros 
of visibility can occur only if the phase 9(V) of AB oscilla- 
tions is independent of bias V. A phase independent of bias 
is ensured by symmetry for an interferometer with L\ = L 2 
and t\ = 1/2. For other parameter choices, our model, in 
which there are no interactions between electrons on opposite 
interferometer arms, yields a bias-dependent 9(V), essentially 
because that Hartree potentials for electrons on each arm vary 
differently with V. For this reason, it would be interesting to 
generalise our treatment to systems with interactions between 
edges as well as within each edge. We expect that such inter- 
actions would reduce or eliminate the dependence of 6(V) on 
V, giving near or exact zeros of visibility even for t\ 7^ 1/2 
or L\ 7^ L 2 , in accord with experiment. Other possible gener- 
alisations include calculation of finite temperature effects and 
of noise power. 
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Appendix A: Calculation of the commutation relations 



We wish to show that the 



commutators 



|«Sm(*i),«*m(*a)] and [G 12 (h),n b tun (t 2 )] appearing 
in Section ITV Bl are zero for t\ > t%. As a first step we obtain 
commutation relations for the fermion fields appearing in 
Eq. (l80t of section HVFl To keep notation concise, we define 

*„ (z v ) = ijj v (d v ,t) = e-if-Mv-^M*)^ {Zv) t 

where z v = d n — vpt. The kernel Q v (z) is nonzero only 
for z > —drj. This follows from causality (with interactions 
occurring only inside the interferometer, the phase of an elec- 
tron at QPC b cannot be influenced by electrons that have yet 
to enter the MZI), and can be demonstrated explicitly for the 
case of charging interactions, using the analytic structure of 
the phase shifts, Eq. d86l >. We have 

(x) , ¥„, (x')} = 0, {4> v (x) , ¥+ (x')} = (Al) 



for any 77, rj and x 1 — x > d v . From this we find 

[v a e ia i>t{-VFh)i>2{-v F h) + h-c, 
v b e ifj iff(di - v F t 2 )y 2 (d 2 - v F t 2 ) + h.c] = 

for t% > t 2 . 

We also require the anticommutators of ^> v (x) with itself 
and with (x). These are straightforward to obtain starting 
from the bosonized form for the fermion operators. However, 
it is instructive also to derive them directly from the represen- 
tation in Eq. d80l l. as follows. We have 

% (x) %> (x') = (x') (x) x 

e -[Qr,,Q v '] e i[Q v (x'-x)-Q^l y x-x')]8, 1 , 1 , 

where 



[Qrjj Qq' 



) v (z - x) Q v > (z' - x') X 

[f>r, (z) , p v > (z')]dzdz'. 



Because of the presence of the filled Fermi sea, the density 
operators do not commute. 35 Instead 



[pr, (z) , p v > (y)j = - -z-dj (z - y) 5r,n> 
Zir 



From this 



[Qr,,Qr,'] = 7^ I dyQ (y ~ x')d y Q (y - x) . 
We now prove that 



Q(x'-x)-Q(x- x') = 



2tt 



+ OC' 



dyQ{y- x') d y Q (y - x) . 



Introducing the Fourier transform of the kernel we have 

Q(x'-x)-Q(x- x') = 

J dqQ{q) [eM*'-*) - eM*-*')] 

and 



dy Q (y - x') d y Q (y - x) 



i [ dqqQ(q)Q(- q )e l < x '- x ). 



2vr 



Thus we wish to show 

Q(q)-Q(-q)=iqQ (q)Q(-q). 
It is easy to check that this is the case, using the explicit form 

Q (q) = 2n 1 df (qd/2)(l + 7 e-^/ 2 jo(< 7 <i/2))- 1 . 
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for the kernel in the case of a charging interaction, where 

jo(x) — (sinx)/x. This gives 

(X')} = 0, 

Or')} - S w 6(x-x'). 
It follows directly that [G 12 (t), H b tun {t')] = for t > t'. 



Appendix B: Calculation of the S matrix 

We require an explicit form for the unitary transformation 
of fermion operators generated by 



5 a (t)=Texp -- / H? un (T)dT 



(Bl) 



The fermion operators appearing in Hf un commute with Hint 
because their position coordinates are before the interact- 
ing region. For this reason the interaction representation of 
Hf Un (t) has the simple form 

K? un (r) = v a e la ^+ (-v F r) fa (-v f t) + h.c. 

The Schrodinger operators tp v (—vpr) anticommute at differ- 
ent values of their argument. As a consequence 

mtun (r),-Ht un (r')} = 0, 

for any r. r' and the time-ordering in ( IBU can be omitted. 
Defining 



% (x,t) = S a+ {t)^ n (x) S a (t). 
and using a Baker-Hausdorff formula 



(B2) 



oo 

Ae B =J2-AAB] n 



= A+[A,B] + -[[A,B],B} + 



we obtain 



tpr, (x,t) 



where 



_ j L, u 5;i (x) , -v F t < x < +0 



tl> v (x) , otherwise 



(B3) 



S a = 



-it a e 



Appendix C: Matrix elements of the exponential operator 

Here we present a derivation of the equation which we use 
to evaluate matrix elements of the form 



C kl = (a\c+e^ M ^ 



(CI) 



where the fermionic operators Cj obey usual anticommutation 
relations {cj, c^} — <5y and {cj, Cj} — 0. The matrix My is 
Hermitian and so the exponential 



(C2) 



is a unitary operator. The matrix elements are calculated 
with respect to the states 



| a) = |mi,m 2 • . . ,mjv), 
= |ni,n 2 , . . .,njv) 



(C3) 
(C4) 



with fermions occupying single particle levels enumerated as 
mi < m2, • ■ • , < to at and 7ii < 712, . . . , < tin correspond- 
ingly, here TV is a total number of electrons. We can write 
these states as a product of creation operators acting on vac- 
uum 



I") = c^ JV c+ JV _ 1 ...c+Juoc) 
I/ 3 ) = ..c+|uac). 



(C5) 
(C6) 



The matrix elements jC 1 b can be written as 

CfeZ = (vac|c mi c m2 . . . c mjv x 

£ fc ^ c+ N c+ N1 ... c+ I vac) (C7) 

or after commuting the operator to the left and cj to the 
right we obtain 



C kl = (-!)*-«■ 



VCLC\C mi C rn2 . . . C rnp _ 1 C rrip j r _ 1 



U ci N ^t w _i ■ ■ ■ _! . . . c+ I vac) (C8) 

where p and g are defined such that k = m p and I — n q . Now 
using the unitarity of the U matrix we can rewrite it as 

Cu = {-l) p+q (vac\c mi c 71 



T+TTf.+ 

TIP, 

S+ 7-f 



U cZ M U + Ucl M M + ...Ucl X1 W 

njv TIN — 1 "■g+i 



E/c+_ l C/+...f/c+<7+t/|va C ) (C9) 
Using the fact that [/ |t? ac) — |i> ac) we arrive at the equation 

Cki — ( — 1 )^ ^ (i'flc| c^7 I , 1 c m2 . . . c rrip _ 1 c mpJrl . . . c m N x 



£ ^c+ JV _ 1 . . . c+ +i c+ _ 1 ... 5+ |«ac) (CIO) 

where we defined = UcllJ + . Applying the Baker- 
Hausdorff identity we obtain 



= J2^ M l M ^--- M L^i = J2 u M ( C11 > 

n=0 ' j=0 

where matrix U is defined as 

U = exp[iM]. (CI 2) 

Substituting Eq. (IC1 U into Eq. ( IClOb gives an equation for 
the matrix elements 

Afc^D^detD (C13) 

where matrix D is formed from the matrix elements with 
indices i and j spanning the occupied states mi, 7712, . . . 
and 77i, 77 2 , . . . 7ijv correspondingly. 
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Appendix D: Numerical evaluation of the correlators 

In this Appendix we provide details of the numerical pro- 
cedure, outlined in Sec. IIV Fl which we use at the final stage 
of our calculations to evaluate the correlation functions from 
Eq. g2). 

We consider a system of length L with periodic bound- 
ary conditions, which leads to momentum quantization k — 
2irn/L, where n E Z. A cutoff is introduced on the num- 
ber of momentum eigenstates, so that n E [—N max , N max \, 
where N max is a positive integer with the total number of 
states N to t — 2N max + 1. The largest value of N max which 
we used in the calculations in order to check the convergence 
was ~ 1000, although in most cases N max ~ 500 was suf- 
ficient. The numerical calculation can be divided into three 
steps. 

Step I. We generate a matrix M with dimensions 2N to t x 
2N tot , where the factor of two originates from the number of 
electron channels in the problem. The structure of this ma- 
trix is as follows. In the matrix element M a p the indices a 
and (3 denote channel and momentum of the creation and an- 
nihilation operators in the kernel. We reserve even indices for 
the first channel and odd for the second channel such that the 
matrix element of the kernel operator between states k and k' 
is represented by 2 x 2 block in the matrix M. The differ- 
ence q = k — k' gives the plasmon momentum q, which is 
substituted into Eq. (l82l with the corresponding form of the 
plasmon phase shifts. These phase shifts are obtained from the 
Eq. (l86l > in the case of charging interactions, Eq. (l93l for ex- 



ponential interactions and can be calculated numerically using 
Bogoliubov equations for other interaction potentials as, for 
example, we have done in Sec IV CI for Coulomb interactions. 
For the zero mode i.e. q = we use Eq. (1531 . 

We evaluate the kernel from Eq. (l82l for every set of k, k' 
for each channel (the kernels are different if arm lengths are 
not equal). Each 2x2 block with given k, k 1 is multiplied 
by the 2 x 2 matrix S^*S^ p or S%*S%p (according to which 
channel it originates from), where S a is defined in the Eq. 
(134-b . The difference of the resulting matrices gives the (k, k') 
block in the matrix M. From this matrix we calculate H = 
exp(iM), which we use in the next step. 

Step II. Next we generate a matrix D with dimensions 
N p x N p , where N p = Np + N% is the total number of 
particles in both channels, and Nff is the number in chan- 
nel [i. Convergence was achieved in most of the cases for 
~ Np ~ 400. In the initial state particles occupy mo- 
mentum eigenstates -^N max . . . ^-(-N max + N^ 2 ). To 
study voltage dependence of the correlators we fix number of 
particles in one channel such that the states {—^N max ... 0) 
are occupied and we vary the number of electrons in another 
channel. The matrix D is obtained by taking matrix elements 
of H which correspond to occupied momentum eigenstates in 
the initial state. 

Step III. We calculate the determinant and the inverse of 
the matrix D. Using this and Eq. ( IC 1 31 > we evaluate corre- 
lators Gkk' in momentum space for every pair of momentum 
indices. Correlators in real space are obtained by performing 
a discrete Fourier transform of Gkk' ■ 
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